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Abstract. A family of maps or flows depending on a parameter v which 
varies in an interval, spans a certain property if along the interval this property 
depends continuously on the parameter and achieves some asymptotic values 
along it. We consider families of periodically forced Hamiltonian systems for 
which the appropriately scaled frequency T3(v) is spanned, namely it covers 
the semi-infinite line [0, oo). Under some natural assumptions on the family 
of flows and its adiabatic limit, we construct a convenient labelling scheme 
for the primary homoclinic orbits which may undergo a countable number of 
bifurcations along this interval. Using this scheme we prove that a properly 
defined flux function is C 1 in v. Combining this proof with previous results 
of RK and Poje, immediately establishes that the flux function and the size 
of the chaotic zone depend on the frequency in a non-monotone fashion for a 
large class of families of Hamiltonian flows. 



1. Introduction 

The modern approach to the study of a dynamical system identifies solving the 
problem with finding the global structure of phase space. On the other hand, bi- 
furcation theory has mostly concentrated on local behavior in the parameter space 
(so even when global phase space bifurcations are considered, they are always con- 
sidered in a small neighborhood of a given parameter value). Here we suggest that 
global analysis with respect to parameters leads to understanding of the phase space 
structure in parameter regimes which are unreachable by the currently known anal- 
ysis techniques. In particular, we examine some features of the global dependence 
of periodically forced systems with finite forcing amplitude on the forcing frequency. 

The study of the effect of periodic forcing on homoclinic loops has been thor- 
oughly investigated since the times of Poincare. Analytical understanding of such 
systems has been attained in three possible limits (see 0,0): the fast and slow 
oscillations limits and the small forcing limit (see figure . For fast oscillations 
averaging methods apply, and together with KAM theory, these are used to obtain 
upper bounds on the separatrix splitting. These bounds are exponentially small 
in the frequency ([17|.[T%]). Lower bounds on the separatrix splitting involve del- 
icate analysis which has been proved, under some structural assumptions on the 
flow, only in recent years (see |S] and references therein). The slow oscillations 
regime corresponds to the region where adiabatic theory applies. Two approaches 
have been applied to this limit - the first one corresponds to extending classical 
adiabatic theory to separatrix crossing (^ELE!]) an d the other corresponds to ex- 
amining the geometrical properties of stable and unstable manifolds to hyperbolic 
manifolds with slowly varying motion on them (see [5] , JU| i [201 i ID ) • Most of the 
analytical studies have been performed in the small oscillations regime in which 
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FIGURE 1. Frequency spanning families. 



perturbation methods apply (see (jHj) and references therein). The most impor- 
tant ingredient from these studies to the current paper is the existence of a flux 
mechanism via lobes (see |2j.[TTj.p ] .p ] .p ] ). 

Here we show that for a large class of systems, the phase space structure and its 
associated transport properties have some common non-monotone dependence on 
the frequency of the forcing term. Namely, we establish that in some cases, even 
for finite size oscillations, there exists a function which depends on the frequency 
continuously. After some natural scaling, this function is simply the sum of the areas 
of the incoming lobes per period, normalized by the forcing period. Hence, results 
obtained in the fast and slow limits supply information regarding the behavior in the 
intermediate regime where no analytical methods apply (see figure^). Furthermore, 
this proposed viewpoint leads to a non-traditional scaling of the flow which is 
relevant for studying its behavior near homoclinic tangles. This paper supplies 
mathematical formulation and generalization of the common work of the author 
with A. Poje |18|. in which similar issues were considered in the context of fluid 
mixing. 

The paper is organized as follows: In the first section we define the class of 
forced Hamiltonian families which we consider. Roughly, these are forced Hamil- 
tonian families with split separatrices, for which the splitting exists for all positive 
parameter values and as the family parameter varies the scaled frequency includes 
the semi-infinite line. In the second section we construct a labelling scheme for 
primary homoclinic orbits, which relies on the concepts of unstable ordering in- 
troduced by Easton 3 . This labelling scheme is valid for all parameter values 
and behaves well across homoclinic bifurcations. In the third section we prove 
the main result of this paper - that the flux through the homoclinic loop depends 
continuously (in fact it is C 1 ) in the spanning parameter of the family (in |18| we 
considered a restricted class of systems for which the continuity of the flux followed 
immediately). To prove this result several properties of the lobes are studied, the 
details are included in appendix 2. In the fourth section we recall the proof of RK 
and Poje, showing that the flux is non-monotone in the scaled frequency and recall 
the implications of these results, especially in view of their implications regarding 
the stochastic zone width (see [22]). Section five includes a demonstration of these 
concepts for the forced center-saddle bifurcation problem. Conclusions are followed 
by two appendices - in the first the homoclinic scaling is discussed whereas in the 
second detailed proofs of some properties of the lobes are included, the relation 
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between the flux and the areas of the exit and entry sets is explained, and the 
adiabatic limit is discussed. 

2. Frequency-spanning families. 

Consider a one- and- a- half degree of freedom Hamiltonian H{x,y,Lo(i/)t + 9;v) 
which is 2n periodic in its last argument and depends smoothly on the parameter 
v € [0,i/*), where v* > may be finite or infinite. This Hamiltonian may be written 
in the form 

(2.1) H(x, y, w{v)t + 9;v) = H (x, y; v) + Hi (x, y, u>{v)t + 9; v) 

where the second term has zero time average: 

Hi (a;, y, oj(i/)t + 9; v)dt = 



and may be as large as the first one. Assume the Hamiltonian system of (|2.1|) : 

dx dH(x,y,u(i/)t + 9) 
^ ' lit = dy = ' V; ' ^ ' y ' ' ' V ' 

dy dH(x,y,u;(i/)t + 9) 

— = — = v Q [x,y;v) +v 1 {x,y,uj{i/)t + 9;i/) 

satisfies the following structural assumptions: 

Al: Hq and Hi are C r ,r > 2 functions of all their arguments, Hi is 2tt 
periodic with zero mean in its last argument . 

A2: For all v values in [0, v*\, (E2> possesses a hyperbolic periodic orbit which 
depends smoothly on v. 

A3: The family of Hamiltonians spans the frequency parameter: the fre- 
quency uj(v) is a smooth C r ,r > 2 function, and it maps the interval 
[0,i/*) onto the semi-infinite line: [0, 00) = {lo(v)\v E [0,i/*)}. Further- 
more, w(v) — > only at one of the interval's boundaries and at no interior 
point. 

A4: For oj = 0, for any phase 9, at least one branch of the stable and un- 
stable manifolds of that periodic orbit coincide to form a homoclinic loop. 
Furthermore, the manifolds and the homoclinic loops created at oj = de- 
pend smoothly on 6, and the phase space area enclosed by these homoclinic 
loops, n(R L ' R (9)), is non-constant and has a finite number of extrema: 

^{R L ' R {9)) ^ 0, for all 9 G [0, 2ir\\{B x , 9 2 , 9 No }. 

For concreteness, with no loss of generality, we assume that at u> = the mani- 
folds topology is of a lying figure eight (closed geometry) or of a fish swimming to 
the left (open geometry), so the left branches always coincide at u> = 0, see figure 

121 

For consider the standard global Poincare map F of the extended phase 

space (x, y, t), which is simply the time-^j- map of the time dependent flow (here 9 
is introduced to account for phase variations): 

27T 2-7T 

F e : {x{Q;0,v),y{0;0,v)) -> (*(— ; 9, v), y(— ; 9, v)) 

The two dimensional symplectic map Fq has, by the above assumptions, a hyper- 
bolic fixed point ^{v) with its associated stable and unstable manifolds, p is a 
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Figure 2. Manifolds geometry at u = 0. 

a. Closed geometry (lying 8) b. Open geometry (fish). 

primary liomoclinic point of 7(f) iff the open segments of the unstable (respec- 
tively stable) manifold connecting with the homoclinic point p, U("f(u),p) 
(respectively S(j(u),p)) do not intersect ( , [3*].|T§]1: 



The points p\ in figure|2|are all primary homoclinic points, whereas the intersection 
points of the manifolds inside the shaded region in figure El are homoclinic points 
which are not primary. We assume: 

A5: For all v ^ there is a finite number of primary homoclinic orbits. 
Furthermore, for each intersecting branches of the manifolds there exist at 
least two primary homoclinic orbits which are topologically transverse. 

A6: The number of primary homoclinic bifurcations is bounded in any bounded 
open interval of v for all v 7^ 0. The order of the primary homoclinic bi- 
furcations is uniformly bounded by the integer M. 



Finally, to establish asymptotic behavior in the frequency, we assume there exists 
an appropriate scaling so that the scaled separatrix length scales and the maximal 
velocity along the separatrix are of order one: 

A7: There exist smooth (C r in v) non- vanishing scaling functions, Wq{v), W\ (u), 
L x {y), L y (y), i max (i^), such that in the scaled variables: 



For convenience of notation, let us assume that M < 10. 



(2.4) 



(2.3) 




H 



H. 



W {v)L x (v)L y {v) 



(2.5) 




the scaled system: 



(2.6) 



dx 

H 



u (x,y,v) + Aui{x,y,wt + 6;v) 



Vo(x,y;y) + Av x {x, y, ut + 9; v) 
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with the two non-dimensional parameters A, u> : 
A(y) = , uo\v) — io(y) 



W {yY w w W (v) ' 

satisfies assumptions A1-A6. Furthermore, for the scaled system, the width 
and length of the separatrix and the maximal velocity along it 1 are all of 
order one for all v £ (0, v*). 

A(v) measures the relative strength of the temporal oscillations to the mean 
flow whereas Z3(y) compares the oscillation's time scale with the travel time along 
the separatrix loop (outside of the saddle orbit neighborhood). We assert that the 
magnitude of these two parameters supplies complete information on the qualita- 
tive behavior of the system (|2.6|l and hence on the original system. Notice that 
assumption A3 is now made with respect to the scaled frequency. 

Definition 1. A family of Hamiltonian systems depending on a parameter v which 
satisfies assumptions A1-A7 is called a frequency-spanning homoclinic family and 
v is called a frequency-spanning parameter. 

A few remarks are now in order: 

• Assumptions A1-A7 clearly hold in the standard near integrable case (small 
A(u)) where the steady flow has a homoclinic loop of fixed size (independent 
of v), the perturbation is generic (so A5 and A6 are satisfied), and its 
frequency uj(v) — v, spans the half real line. 

• A5 refers only to primary homoclinic bifurcations - otherwise it would have 
been violated generically, see and references therein. 

• In the ncar-integrable case the scaling functions may be simply extracted 
from the integrable flow, see for example section [S] In appendix 1, after 
some notation is established in the next section, we propose an algorithm 
for defining the characteristic scales in the finite amplitude size forcing case. 

• The assumptions here are slightly weaker than the ones in |18| : in particu- 
lar A5 and A6 replace the stronger assumption of JHj that there exists one 
topologically transverse primary homoclinic orbit which depends continu- 
ously on v, a property which may be easily violated in applications. Relax- 
ing this assumption requires the introduction of a new labelling scheme for 
the homoclinic points which takes into account the possible annihilation of 
primary homoclinic points. 

Hereafter we assume the forced system is in its scaled form and we drop all the 
over bars when non- ambiguous, as described next. 



3. LABELLING SCHEME FOR THE PRIMARY HOMOCLINIC POINTS 

Denote by P{v) = {p r " {lJ \v),n = 1, N(f)} the topologically transverse pri- 
mary homoclinic orbits of 7(f) belonging to the left branch (for simplicity of no- 
tation we consider hereafter only the left branch of the loop), where p Tn ^>[y) = 

{Pi, n (^)}iS:-oo = {-f ,t Po"^( I/ )}^-c») namely the discrete index i denotes iter- 
ations under the map and the index r n (v) denotes the label of the specific orbit. 
Each such transverse homoclinic point, p r ™(y), defines a region, R^ n (v), enclosed 



See appendix 1 for the precise definitions, which use some of the notation of the next section. 



6 



VERED ROM-KEDAR. 



by the closed segments of the stable and unstable manifolds which emanate from 
the fixed point and meet at p[" (y) , namely 

(3.1) dI$»(v) = UMv),p?{u)]USky(v),p:»(u)], 
where, clearly 

Denote by < V\ < v% < .... the ordered parameter values at which primary 
homoclinic bifurcations occur (by A6 this may be done), and by T the set of all 
primary homoclinic bifurcation values in [0, v*)\ T = { v^, ■■■}• In figures [3] and 
0] we illustrate a possible behavior of the manifolds before and after a primary ho- 
moclinic bifurcation occurs. We say that p rn ^"> undergoes a homoclinic bifurcation 
at Vi if p r "( Vl ) is a tangent homoclinic orbit. Next we construct a labelling of the 
orbits in P{v) which is convenient even across homoclinic bifurcations. 

Lemma 1. Consider a frequency spanning homoclinic family. Then, there exists a 
labelling scheme for the primary homoclinic orbits satisfying the following proper- 
ties: 

(1) r(y) is piecewise constant, changing only at values of v at which p r ^(y) 
undergoes a homoclinic bifurcation. 

(2) The labelling scheme respects the unstable ordering, so that r n (y) < r m {y) 

fPo (") <uPo ( v )> namely (see 

(3.2) U[-y{v),pp(v)] CC/[ 7 M,pS"»] for r n < r m . 

Proof. By construction. We construct such a labelling scheme inductively, in the 
intervals (yu, Vk+i), starting with the interval (vq — 0, Vx). 

By A5, N(v) > 2 for all v 5^ 0. Since the Poincare. map is orientation preserving 
and since only topologically transverse intersections are counted, N{v) is even. By 
assumption A6, there exists a well defined limit of the number of primary homoclinic 
orbits as v — > + , so that N(0 + ) is well defined and is finite, and v\ > 0. For 
< v < vx, let r n (v) = n,n = 1, iV(0 + ). The labelling is chosen to obey the 
unstable ordering Pq(v) < u Poi v ) f° r i < h where we identify: 

(3.3) p r NM+1 (1/) = pY (i/) - Fp r > (y), and p r ° (v) = p r _\^ („) = F^p^ <y). 

These rules determine the labelling for < v < v\ up to a cyclic permutation, 
and do not define the origin (which homoclinic point is p\(v)\ To remove the 
first ambiguity, define the labelling so that as v — > + , the first homoclinic orbit 
defines the region of maximal area 2 : /i(i?o(0 + )) = max„(/i(i?Q (0 + )). To remove 
the second ambiguity, we define the origin (i = 0) as the point for which the 
boundaries of R}(0 + ) are of minimal length; Denote the arc length of the boundary 
of by \dR^ n (v)\. Notice that |9i?["(^)| — > 00 as i — > ±00. Let Pq(v), satisfy 

\dRQ(v)\ = mini 1^-^1(^)1- When there are no degeneracies (here, when the maximal 
area is achieved at a unique orbit and the minimal boundary length is achieved at 
a unique homoclinic point) this procedure determines uniquely the labelling of the 
topologically transverse homoclinic orbits for < v < v\. If a degeneracy occurs, 
choose any one of the finite number of the maximizing (respectively minimizing) 
orbits (respectively points). 

Given the labelling Lbl{v^,_^) = {ri,r 2 , rjv(„)}, for k > 1, so that Lbl{v^_ x ) = 
Lbl{v) for v S i/fc), we construct Lbl{v^) as follows: 

2 By area preservation ^(R^ n (y)) = ^(-Rq™ (v)) for all i. 
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Figure 3. Primary homoclinic points. 
Here N(v) = 2, v < v\ of figure [SJ 

Case 1. I] the associated homoclinic orbit, Pi"(f), depends continuously on v at 
Vk leave r n unchanged. 

Case 2. Ifp^ n (y) cease to exist for v > Vk , delete r n from Lbl(ut). 

Case 3. If M additional homoclinic orbits are created between (along the unstable 
ordering) p T ™ (y) andp r ^ r>+1 (y), insert the M labels, r J = r n +j-10~ , j — 1, M , so 
that Lbl{v+) = {n, ra, .., r», r„ + Kr fe , r„ + 2- 10" fc , r n + M • 10~ k , r n+1 , r N(v) } 

By construction the labelling scheme respects the unstable ordering and does 
not change labels unless homoclinic bifurcation occur, as claimed. Notice that the 
first case in the proof includes two possibility which are topologically equivalent 
- the first is that the homoclinic orbit p^iyk) is not involved in the homoclinic 
bifurcation and the second possibility is that p^iyk) is a tangent periodic orbit 
with odd order tangency so that no topological bifurcation occurs at v^. □ 

For any v we have constructed a finite, ordered set of labels Lbl(v) — {ri, r-2, r N ^} 
(so ri < r 2 < .. < r N r u )) with corresponding homoclinic orbits which satisfy H3.2[l . 
The ordering k{r{v)),r £ Lbl(v) enumerates the primary homoclinic orbits in P(y) 
by their unstable ordering. While the r' s depend smoothly on v as long as p\ (y) has 
not bifurcated, k(r{v)) can be discontinuous in r{v) due to bifurcations of orbits 
p[ (v) with some r' < r. 

For example, in the simplest, generic case of quadratic tangency, if at V\ two pri- 
mary homoclinic orbits are created between (in the sense of the unstable ordering) 
the jth and the j + 1 homoclinic orbits which existed for v < v\ , then these will be 
labeled asr= = j.2, so that p 3 Q < u p^' 1 < u p^' 2 < u p^ +1 (^), as demonstrated 

in figure 0] 
(3.4) 

P{y)U<u<* = {pH"),p 2 ("), ...yHy u H,/ 2 My +1 w, ..., P N (° + \v)}. 

By construction, for < v < v\, = k, or equivalently, k(r;v) — r. On the 
other hand, for v\ < v < v 2 , for the generic case of l|3.4[l we have Lbl(v)\ Vl<v<V2 = 
{l,2,..,j,j.l,j.2,j + l,...,N{0+)} = {r 1 ,..,r N{0+)+2 }so 

!r for r < j 

[r] + [(r — [r]) • 10] for j < r < j + 1 for v\ < v < v% 
r + 2 for r > j + 1, 
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Figure 4. labelling after the first bifurcation. 
v =v_ where v\ <v< V2 sec figure [5] 




VI V V2 V3 V4 



Figure 5. Primary homoclinic bifurcation diagram (schematic). 
t represents natural parametarization along the unstable manifold. 

and we see that k(r; v) changes at v\ for all k > j, hence, as opposed to the 
r„'s, it does not serve as a good labelling system for the homoclinic orbits across 
bifurcations. 

4. Lobes topology 

Consider now a fixed v value, with Lbl{y) = {r\,r2, ■■■,r2n} where N(u) = 2n, 
and Pq 2 " +1 — FpQ 1 = p^ 1 ■ Henceforth we omit the explicit dependence on v unless 
needed. 

Definition 2. p r ,p r G P{v) are called neighbors if \k(r) — k{r')\ = 1. 

Definition 3. Rq and Rq (see are called neighboring regions if p r ,p r are 

neighbors. 

The difference between neighboring regions i?Q and Rq are lobes - regions which 
are enclosed by the segments U[pq,Pq ]U S[pq,Pq ]: 

Definition 4. A lobe L{p r il p\ ) is the region enclosed by dL = U\p r il p\ ]U S\p r i ,p\ ] 
where p\-,p\ are neighboring primary homoclinic points. 

4.1. Incoming and outgoing lobes. Let 

Dk = L{p T Q k ~ 2 ,p T Q k ~ 1 ), E k =L(p r ^-\ P ^), fi rfc=l,...,JV(i/)/2 
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where (lrt.3|) is used to define D%, see figure 01 Next we prove that this definition 
indeed corresponds to the figure, namely that the choice of the labelling is such 
that -Dfc+i is exterior to i?Q 2fc whereas Ek is inside the region i?Q 2fe_1 : 

Lemma 2. For all v > and all k the lobes Ek and Dk satisfy 3 : 

(4.1) R^ 2 "- 1 = R r ™ M E k , Rl 2 "' 1 = Rq 2 "' 2 W -Dfc . 

Proof. First we prove that for < v < i>\, the lemma follows from the definition 
of the labelling and claim ^ which is proved in appendix 2; recall that for these v's 
rk = k and Rq has maximal area. By the lobes and regions definitions, it follows 
immediately that neighboring regions defer from each other by lobes, so proving 
the lemma amounts to proving that: 

(1) The lobes boundaries cannot intersect the corresponding boundaries of the 
regions, hence they can be either completely interior or completely exterior 
to the regions. 

(2) The definition of the exiting (D's) and entering (E's) lobes is consistent 

with jnj. 

The first part is contained in claim ^ of appendix 2, where it is proved that: 

2# +1 )rh T dB% =0, 

(4.2) S(p? +1 ,^)rihT^ +1 =0, 

where ftly denotes topologically transverse intersection. These results, together with 
the observation that S\pq{v),Pq(v)\ C 8Rq(v), implies that the lobe E\, which is 
enclosed by the segments U[pq(v),p^(v)\J S\pq{v),Pq{v)\ is either contained in Rq 
or is completely outside of it. Hence, either Rq = Rq W E\ or Rq = i?J ttl E\. 
However, the latter contradicts the choice of R\ as the region with maximal area 
as v — > + . Another manifestation of the labelling scheme may be formulated by 
looking at the manifold's orientation: E\ C R\ iff U\jPq(v),Pq(v)\ is interior to 
S[pq(v) , Pq(v)] as shown in figure|3l The manifolds orientation is clearly preserved 
for v £ (0, v\). Furthermore, by topological transversality of the points in P(y) it 
follows immediately that all odd (respectively even) indexed points have the same 
manifolds orientation as of p\{v) (respectively Pq(v)), so that for < v < v\ the 
lemma is proved for all k. 

For Vj <v < Vj+i, the lemma is proved by induction, noticing that at the bifur- 
cations, topologically transverse homoclinic points are inserted/deleted in neigh- 
boring pairs, hence the orientation of the manifolds at PQ 2k+1 (v) (respectively at 
p T Q 2k {v)) is always as that of Pq{0 + ) (respectively as of p§(0 + )). Notice that tan- 
gent, non topologically transverse homoclinic points are not labeled. □ 

Corollary 1. For all v > the regions with even order have smaller areas than 
their neighbors, namely i?Q 2fc+1 are regions of locally maximal area whereas R r Q 2k are 
regions of locally minimal areas. 

Finally, note the following property of the lobes: 

Proposition 1. The lobes Ek,Dk k = 1, ...,N(y) are simply connected. 



'I+J denotes a union of sets with disjoint interior. 
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Proof. The lobe's boundaries, by definition, are given by U{p r 3 (v),p 3+1 {v)) U 
S(Po ( v \Po +1 which are both connected arcs which are joined at their end- 
points, hence the lobes are connected. Each component of the boundary cannot 
intersect itself, and by lemma^they cannot intersect transversely each other either, 
so the only source of non-trivial topology may stem from tangencies. The appear- 
ance of primary homoclinic tangency (e.g. at v = v\ of figure^) does not change 
the simple connectedness property of a lobe. Finally, it follows from eq. (|4.2|l that 
the lobe cannot reconnect and form a ring (even in the closed geometry case!). □ 

4.2. Continuous dependence on parameters. For every v there exist a collec- 
tion of incoming and outgoing lobes. Their areas changes with v, and, at homoclinic 
bifurcation points these changes can be discontinuous (at even order bifurcations 
where lobes split or coincide, as shown in figure EJ. Nonetheless, 

Proposition 2. Consider a frequency- spanning homoclinic family. The sums of 
the areas of the outgoing and incoming lobes: 



are C 1 in v. 

Proof. It follows from lemma ^ that away of the primary homoclinic bifurcation 
values {vi]°^ 1 the area of each lobe depends smoothly (C r ) on v, hence so do 
fj,D and fiE- By orientation preservation, near the bifurcation values an even 
number of topological transverse primary homoclinic points are added or elimi- 
nated. Consider first a tangential bifurcation of even order occurring at Vk so that 
2n new topologically transverse homoclinic points are created at v = . These 
new homoclinic points, appearing, for example, between the homoclinic points 
p 2 - 5 ,p 2 ^ +1 , are denoted by p 2j m , m = 1, 2n see figure 5 4 . The bifurcation splits 
the lobe D } {v x ) (enclosed by U[p 2j ,p 2j+1 ] U S[p 2j+1 ,p 2j ]) to the 2n + 1 disjoint 5 
lobes Dj(u), Ej + \ (v), -Dj+i (v), Ej +n {v), Dj +n (u) created by the new homoclinic 
points. As v — ► v^, by the smooth dependence of the manifolds on v (recall that 
only a finite extension of the local stable and unstable manifolds is considered here) , 
the areas of the interior, newly created lobes (£j+i(^), Ej+ n (v)) vanish, whereas 
n(Dj(v)) + ... + n{D.j +n {v)) -^ v ^ v + /i(-Dj(fi)). The area of the diminishing lobes 

near is of the form (y — Vk) N+1 ^ N where N is the order of the tangency at 
Vk- Clearly the same argument applies to the case where pairs of homoclinic orbits 
annihilate at v^. It follows that generically ^i(D) and n(E) are C 1 in v (I thank B. 
Fiedler for this observation). □ 

Intuitively, we think of the lobes as existing lobes and the Ek as entering 
lobes. A more elegant framework for discussing transport in phase space is achieved 
by introducing the notion of exit and entry sets (see [lSp: 



^The other case, of homoclinic bifurcation occuring between the homoclinic points p 2 -' - 1 ,p 2j , 
is similar: simply change the corresponding indices and interchange the letters E and D in this 
paragraph. 

5 by smoothness of the manifolds 
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Definition 5. For any region R, the exit (Dr) and entry (Er) sets of R are defined 
as: 



In appendix 2 we study the relation between entry and exit sets of the regions 
Rq 3 and the corresponding entering and exiting lobes. In particular, we define, for 
each region Rq 3 the set of turnstile lobes , in terms of union of entering and 
exiting lobes (see equations l|9.6[) in the appendix) . We conclude that in case there 
are no intersections between these two sets they correspond exactly to the entry 
and exit sets of Rq 3 . However, if such intersections exist, the two sets are different. 
We provide an example which demonstrates that when such intersections exist the 
measure of the entry and exit sets depends on j. 

We therefore conclude that while the notion of entry and exit sets appears more 
elegant, it has a major flaw from a dynamic point of view - it is not invariant with 
respect to the regions definition. 

Another problem arising from the definition of the entry and exit sets has to 
do with self intersection of the lobes (for close geometry, in the limit of highly 
stretched lobes). Indeed, the boundaries of the entry and exit sets may depend on 
the appearance of non- primary homoclinic orbits. The appearance of such orbits is 
expected to depend sensitively on parameter values. Equivalently, the proposition 
regarding the smooth dependence of the lobes area on parameters is proved for the 
sum of their areas, \ie I'see l4.3|) . and not for the area of the lobes union. In fact, any 
definition of flux which depends on homoclinic orbits of higher order (namely not 
primary homoclinic orbits) must be carefully examined due to the recent results on 
Richness of chaos [23] in the neighborhood of a homoclinic tangency. 

Hence, we define the flux into a region using the concept of the sum of lobes 
areas: 

Definition 6. Consider a frequency spanning homoclinic family. The scaled flux 
function of this family is: 



Corollary 2. The flux function of a frequency spanning homoclinic family depends 
continuously and has a continuous derivative (is C 1 ) in the frequency spanning 
parameter. 



Notice that the scaled system defines a scaled flux function. Substituting the 
scaling function, we observe that the flux function is scaled by the typical area and 
time scale associated with travelling along the homoclinic loop, as appropriate: 



D R = R\F~ 1 R, E R = R\FR. 



(4.4) 




Proof. By proposition [3 and assumptions A3 and A7. 



□ 



flux(v) 



2tt 




Wq(v) L x (y)L y {v) Areai, 



oop 



flux{v) 



5. NON-MONOTONICITY OF THE FLUX FUNCTION 



Theorem 1. Consider a frequency- spanning homoclinic family. Then, in the scaled 
system \2. 6)) the scaled flux function ( eq ) depends non-monotonically on the 
scaled frequency zu(v). 
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Proof. Recall that by corollary [3 the flux function depends continuously on v and 
that fiux{v) = ^rT^E- As v — > 0, the system (|2.6|) approaches the adiabatic limit. 
Then, by corollary |3| of Appendix 2, ^e{ v ) — * Me(0 + ) > 0, which implies, by 
assumption A3 and A7, that j jiuxiy) — ► 0. Furthermore, since V\ (the hrst primary 
homoclinic bifurcation value) is bounded away from zero, it follows from A3 that 
there exists a < v < v\ such that ffi U x(y) ^ and ffiuxty) is monotone in 
the interval (0, v\. On the other hand, as w(y) — * oo, averaging of the scaled 
system Ij2.fil) implies that the separatrix splitting is exponentially small in \jw(y) 
(see It follows that he{v) — ► exponentially as v — ► v°° (where 

linv-^oo zu{v) — oo, and by A3 v°° 6 (0, f*)), and so does jfi U x(y)- In particular, 
there exists a. v' < v* such that for > v', ffi ux (v) < ffiuxip)- This proves that 



ffiuxiy) is non-monotone in v. □ 

We now discuss some of the implications of the non-monotone behavior of the flux 
function, see JS| for a fluid mixing application of these results. First, we notice that 
in the near-integrable case, when the number of primary homoclinic orbits is fixed, 
the flux function is proportional to the amplitude of the Melnikov integral. Hence 
the above theorem proves that the Melnikov function amplitude is non-monotone 
in this case, which indeed explains the non-monotone figure one typically gets in 
numerous calculations of the Melnikov integral (e.g. forced Duffing's, forced Cubic 
potential, forced pendulum). It is clear that very small (large) flux corresponds to 
a small (large) chaotic region near the separatrix, hence, it was suggested that the 
Melnikov function amplitude gives a good characteristic to the amount of chaos 
in the system. However, we proved that the flux function is non-monotone in the 
spanning parameter. It is therefore natural to compare the dynamics, and the 
properties of the chaotic region at equi-flux parameter values. Let v a ^ Uf, denote 
cqui-flux parameter values such that: 



Iiux{v a ) = fi ux {v b ) and w{v a ) < m(v b ) 
It follows immediately, by the flux definition, that 



namely, the sum of the lobes areas of the lower frequency is larger than that of 
the higher frequency. This observation, which is well reflected when one examines 
Poincare maps of equi-flux frequencies may lead to the wrong conclusion that adi- 
abatic mixing leads to a larger chaotic zone. In fact, at least for small amplitudes, 
we can use [221 to show that the converse result hold: 

l^stoch 

where JJ^toch denotes the area of the mixing zone (the area enclosed by the KAM 
tori which are closest to the separatrix on either side of the loops) . A more precise 
statement of these observations, with numerical demonstrations, may be found in 
|18j . The seemingly contradicting statements are well understood when one realizes 
that for small w the lobes overlap considerably - in the adiabatic limit, the measure 
of overlap between the turnstile lobes may approach the full lobe's area (see [T%] 
for a proof), namely there exists a j such that 

where E J , D 3 are defined by eq. (|9.6(l in appendix 2. 
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6. Forced Saddle-center bifurcation: 
As an example consider the forced saddle center bifurcation: 

x = p 

p = v — x 2 + bx n sinwi. 

For < v <C 1. Near x = —\pv a hyperbolic periodic orbit appears. Its separa- 
trix length scales (width and height) in the (x,p) directions are {L x (y),L p (v)) — 
(0(y/V), 0(z^ 3 / 4 )), so L max = L x . The maximal velocity along the loop scales like 
W = U = 0(z/ 3 / 4 ). Finally, since Wi = Vx^f^ = 0(-^=). we obtain the scaling: 

Hence, the effective homoclinic forcing amplitude and frequency are: 
6- 1 A = ~i — ■n>&=—= 



and the scaled system (with over bars dropped where non ambiguous) is: 

dx 

W Tt= p 

= 1 - x 2 + Ax n sinUi. 

dt 

Namely, the effective frequency and damping increases inversely with the bifurcation 
parameter. For fixed u> and fixed A we recover the well known result that the 
separatrix splitting is exponentially small as the bifurcation is approached. For 
fixed b we see that the exact dependence on n matters. For large 57, an appropriate 
(see [5]) Melnikov calculation for (|6.2|) will lead to a function of the form: 

„ ,, . IjJ . b TTU! Ult 

M(t; is, lo, 6, n) = C{^=)-^^ ex P(~^) «» 

where C(-) is, up to exponential small corrections, a polynomial. If n = and 
b = Pu, then C(^) = K^ (see 0). 

Let b — /3^ 1_n / 2 , and consider a frequency ui(v) which vanishes at some finite 
distance from the bifurcation value and is non decreasing as the bifurcation value 
is approached. For example, take ui = (e — v) a v^ where j3 < 1, e < l,a > 
0, —a < fi < 0. Then, the family of forced saddle-center bifurcation with v S (0, e] 
is a spanning family, satisfying assumptions A1-A4. Verifying assumptions A5- 
A6 amounts to obtaining a lower bound on the separatrix splitting for equation 
(16.21) as uj — > oo, a task which has been achieved in [S] for the case n = 0, and 
we will assume here that a similar lower bound may be found for n > so that 
A5 and A6 are satisfied. A7 is clearly satisfied with the scaled parameters 16.11 
Furthermore, the scaled frequency uj = (e — v) a v^ x l^ is monotone for v € (0,e], 
and 57((0, e]) = [0, oo) . In particular, taking n — 2, a = = shows that the 
system: 

(6.3) x = p 

p = v — x 2 + bx 2 sin(e — v)t 
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is transformed to the scaled system (16.2(1 with A = b and oj(v) = i^ 1 ■ We 
therefore conclude that the homoclinic structure of this family behaves in a non- 
monotone fashion - there exists at least one value of v, v m € (0, e) such that the 
scaled homoclinic flux function: 

hvJy) = , n ° op fiuxty) = v~^fiux(y) = E / y^K E i{ v )) 

Areaioop 2n V 2 

has a local extrema at v rn . In particular, this implies that the corresponding Mel- 
nikov function amplitude is non-monotone in v. One expects that uj(v m ) — 0(1), 
hence v rn as e 4 , namely the maxima location scales with the size of the interval on 
which the family is spanning. 

7. Discussion: 

We have demonstrated that the concept of spanning families, in which one con- 
siders global dependence on parameters, leads to non-trivial results regarding the 
properties of forced systems. This approach leads to the development of a new tool 
- a novel way of labelling of primary homoclinic points across bifurcations. Using 
this tool we are able to prove several results regarding the global dependence of 
the lobes on the spanning parameter. We proved that lobes are simply connected. 
We proved that under natural assumptions on the spanning families the flux is 
continuous and non-monotone in the spanning parameter. We demonstrated that 
under the same conditions the area of the exit and entry sets may depend on the 
spanning parameter discontinuously. We demonstrated that this approach leads to 
non-standard view of the forced saddle-center bifurcation. 

The relation of this work to higher dimensional systems, in which the spanning 
parameter is in fact a slowly varying variable is yet to be explored. 

Acknowledgement 1. Discussions with B. Fiedler and D. Turaev are greatly 
appreciated. 

8. Appendix 1: Homoclinic scaling 

Below we supply an algorithm for determining the length and time scales charac- 
teristics for the finite amplitude forcing case. Consider the original system. Define: 

(l x (i,n;v),L(i,n;v)) = I max \x - x'\, max \y-y'\) 

where the regions (v) are the regions defined by the primary homoclinic points 
p^"(v) as in section The characteristic width and length of the separatrix loop 
are given by: 

(L x (v),L y (v)) = [l x {i m ,n m \v),l v {i m ,n m ;v)), 

where 

\\(la;(im,n m ;v),ly(im,n m ;v))\\ = min \\(l x (i,n;v),l v (i,n;v))\\ 

\i\<K,n— l,...,N(v) 

and K is some finite number. Let 

(v) = ma,x{L x (i/),L y (v)} 

The choice of the labelling so that |9i?Q| is minimal among all \dR}\ and the oscilla- 
tory nature of the boundary for large i, implies that i m is expected to be small, and 



FREQUENCY SPANNING HOMOCLINIC FAMILIES. 



therefore independent of K if K is sufficiently large. The characteristic velocities 
along the separatrix loops are similarly defined: 



[Ui{v), Vi(i/)) = limsup |ui(-)|,limsup \vi(-)\ 

v ' \ n,6,t n R t 

Wq(v) = L max {v) max 



n,6,t n,6,t 

L x (v) ' L y {v) 

ih(y) Vx{v) 



,i=l,2 



Wi{v) = L max (p) max ■ 

yL x (y) Ly{v) J 

where Pq" (.) denotes the primary homoclinic orbit with the initial condition p(to \ to) = 
p r Q "(.). Since pjj are all primary homoclinic orbits (in particular, bounded orbits) 
the functions Ui(v), Vi(v), Wi{v) are clearly well defined and are independent of the 
choice of the labelling. Away from the homoclinic bifurcation points L Xl y(is) de- 
pend smoothly on v and i ma x( I/ )) Ui(y), Vi(v), Wi{v) depend continuously on v. At 
Vj these functions may be discontinuous. Let L xy (y) 1 L(y),Vi{y),V i {y),'Wi(v) de- 
note C r - functions, which are close in the C° topology to the corresponding hatted 
functions, uniformly in v. These are the proposed scaling functions for assumption 
A7. 

9. Appendix 2: Lobes 

9.1. Neighboring regions differ by lobes. Here we prove claimQlwhich is used 
to prove lemma |21 

Claim 1. For all neighboring homoclinic points p^ 3 ,p r Q 3+1 € P(y), either U (pq 3 ,Pq 3+1 ) C 
R r Q 3 or U(p r 3 ,p r 3+1 ) n R r Q 3 = 0, namely 

(9.1) U{pl\p r 3+1 )fo T dRl 3 =% 
and similarly 

(9.2) S{p r 3 +\p r J)h T dRl 3+1 =% 

(where A\t denotes topologically transverse intersection). In particular, 

(9.3) U(pl i ,p^)^ T S(^ + \p^) = t 

Proof. Recall that by definition ORq 3 = U^^p^ 3 ] U S[~f,p r Q 3 ]. Since the unstable 
segments cannot intersect each other clearly U(pq 3 ,p 1 Q 3+1 ) n [/[7,Pq j ] = 0. So, to 
prove Q9.1II we need to prove that U(pq 3 ,p 1 3+1 ) rhr SITjPo 3 ] = Notice that 
^[liPo 3 ] = ^[1tPq +1 ] ^ S[pq 3+1 ,Pq 3 ]. Since Pq{v) are primary homoclinic points 
J7(7,Po) nS 1 (7, Pq) = and in particular 

(9.4) u(p r 3 ,p r 3+1 )ns( 1 , P r Q 3+1 )=$, 

and 

(9.5) U( 1 ,p r Q 3 )nS( llP r Q 3 ) =0 

Hence, it is left to prove (|9.3|l . Assume the contrary: denote by Q all the topologi- 
cally transverse homoclinic points in this intersection. Denote by q € Q the homo- 
clinic point which is closest to p 3 by the unstable ordering, namely U(p 3 ,q) i+It 
S(Po +1 ,Pq ) = 0- Then, it follows from (|9.5|l and the observation that S(pq +1 , q) C 
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S(p 3+1 ,Pq) C 5(7, Pq j ) that U(j, q) rfir S(<y,q) = 0, namely q is a topologi- 
cally transverse primary intersection point, which contradicts the statement that 
Pq j ,Pq' +1 are neighbors in P(v). Hence (|9.1|l and i|9.3[l are proven. To prove 
JH3 notice that dR r 3+1 = t/[7,pg J+1 ] U S[y,p T Q i+1 }, hence, to prove the claim 
we need to prove that S(p r j+1 ,p r j ) ftl T U^^^ 1 ] = S(p r Q j+1 ,p r j ) rh T (U[~/,p r Q 3 ] U 
U(pq' ,Pq j+1 )) = 0- The intersection with the first term is empty by (|9.5|l and the 
transverse intersection with the second term is empty by (|9.3|) . □ 



Corollary 3. U(& n i# = iff tf(p? + >),*# + >)) C /or all 

j = l,...,N(u). 

Proof. This follows from the topological transversality at Pq j+1 {y) and the previous 
lemma. □ 



9.2. Exit and entry sets. For any region R, the exit (D R ) and entry (Er) sets 
of i? are defined as (see |12j): 

D R = RXF-iR, E R = R\FR. 

where over bar denotes closure. 

Consider first the case N(v) = 2, where only two topologically transverse primary 
homoclinic orbits exist. From (|4.1|) we immediately get that: 

Rl = Rl W Ex = F^Rl W D x , 

Rq = Rq\ e i 
F 1 Rq = Rq\Di 

hence 

E R x = Rl\FRl = E 1 \FD 1 = Ei\(FD± n Ex) 

D R i = F-\FRl\R^) = F-\FDx\Ex) = Dx\{D x n F~ 1 E 1 ) 

and 

^2 = F(F~ 1 Rl\Rl) = F(Ex\Dx) = FEx\F(Dx n £ x ) 
D Hg = fl^F- 1 ^ = = £>i\(£>i n Ei) 

therefor, we see that D R i = D R 2 = D\ , E R i = Ex and E R 2 = FEx iff Dx fl E x = 
FDx n Ex =0- In particular, in figure HO we demonstrate that if fi(Ex flfli) > 
and Ex fl F.D1 = we obtain: 

fi(E Rh ) = M (i? 1 ) > M (^) = /x^r) - M (£>! n £?!). 

so the measure of the entry and exit sets depend on the region one chooses, while 
the measure of the entry and exit lobes are independent of this arbitrary choice. 
Furthermore, we see that the entry and exit sets depend on the nature of non- 
primary homoclinic orbits. The appearance of such orbits may depend sensitively 
on parameter values. 




Figure 6. Exit sets and lobe. 



For a general value of N{v) (recall that N{v) is even) we define the entering and 
exiting lobes as: 











(9.6) 


D 2k 


= £) 2k ~ 1 


= U D,= U D i U 
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^1+k 
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E 2k 
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E i = U £ i u U f£ i 
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^ fc-1 


(9.8) 
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^ = U '••'.< U >- i: - 











3=1 



It follows (after some manipulation of H4.lt 1 that in the open geometry case 
(where Ei n Ej = Di n Dj = for all i ^ j) the exit and entry sets of the regions 
i?o J are given by: 

E r } =Ei\(Ei riD*) 
= F- 1 (Di\ (E* n Di)) if ^ n ^ = A ° ^ : = for a11 * ^ * 

In particular, in the near-integrable case, where -E J n D 3 = for all j = 1, iV(z/) 
the exit and entry sets are given by the entering and exiting lobes. Indeed, for 
N(v) = 2 we find: 

D 2 = D 1 =FD ll E 1 =E ll E 2 ^FE u 
and therefore the first example follows this general rule. 



9.3. Lobes in the adiabatic limit. The limit of TD(v) — > is called the adi- 
abatic limit. It corresponds to a forcing with frequency which is much smaller 
than the time scale associated with the flow along the separatrix. The tradi- 
tional way of analyzing this limit is by using the adiabatic approximation. In this 
approximation one considers solutions to the frozen Hamiltonian H (x, y, 9; 0) = 
Hq(x,d;§) + Hi(x,y,9;0), namely the time-dependent oscillations are now fixed 
at the frozen phase 6. One then studies adiabatic variables (to leading order, the 
action) , and asserts that in the nonlinear long as the solutions are bounded 

away from the separatrix, such variables cannot change by much due to the per- 
sistence of KAM tori. Furthermore, when the solution crosses the separatrix one 
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can estimate the resulting change in the adiabatic invariances and the phases by 
constructing an adiabatic separatrix mapping (see for example ^H] > [21] > [Hj , [Z] ) ■ An- 
other point of view which is applicable to the same limit (see 0,^01 an d 1201 )j is 
concerned with the geometry of the manifolds and lobes. Consider the extended 
system (x, y, r = uj(v)t), and realize that the manifold {(x, y, t)\x = Xf,y = yf, t € 
[0, 2tt)} is a normally hyperbolic manifold for uj(v) — 0, hence it persists with 
its stable and unstable manifolds for sufficiently small oj(v). The solutions be- 
longing to these manifolds can be well approximated, for semi-infinite time in- 
tervals by the extended system at uj{v) = 0. Denoting the homoclinic solution 
to this extended frozen system by qh(t;To) = (xh(t;To),yh(t;To),To) and the or- 
bits belonging to the stable and unstable manifolds of the extended system by 
q s ' u (t; t , v) — (x s ' u (t; t , v), y s ' u (t; r , v), r = uO{v)t + T ) we see that as W{y) — > + 

q s {t;T 0l is) = q h (t;T Q + uj(v)t) + t G (-co, T] 

q u (t;T ,v) = q h (t;T +ZJ(v)t) + te [-T,oo) 

where dots stand for higher order terms in uJ{y), T is a positive constant and the 
higher order terms are small for W{v) < oj* (T) . Recall that in the limit v — > + 
(which, by A3 corresponds to the limit T3{y) — > + ), there are N(0 + ) topologically 
transverse primary homoclinic orbits, p n = p n (t;TQ,v),n = 1, ...,N(0 + ). By def- 
inition, these orbits belong to both q s (t]T§,v) and q u (t;TQ ,v) for some phase r, 
so that g s (0;TQ,z/) = q u (0; Tq , v). We conclude that these phases, Tq , are given, 
to leading order in uj{y), by the phases at which the area enclosed by the frozen 
homoclinic loop achieves its maxima or minima: 

Lemma 3. (Kaper and Kovacic) Let {9i}^i ' denote the 9's at the minima and 
maxima 6 of the homoclinic loop area /j,(R(9)), where /j,(R(9 1 )) = maxg fi(R(9)), 
and N(0 + ) > 2. Then, as uj(v) — > + , the primary homoclinic orbits p n (t; Tq,u) — > 
Qh(t', On + uj(v)t) and: 

^R(9 n )) = m(-Ro) + o(l) 
where Rq denotes the region enclosed by U[y,p n ] and S[y,p n ]. 

Proof. See where it is proved that to leading order the adiabatic Melnikov 
function, which measures the distance between the stable and unstable manifolds 
is given, to leading order, by 4s(j,(R(6)). □ 

Corollary 4. As oj(v) — » 0, 

li{E k ) = n(R(6 2 k-i)) - ti{R{02k)) 

(9.9) 

hence, provided is satisfied 



The structure of the lobes in this limit is more complicated than the commonly 
seen graphs of homoclinic tangles. In particular, on one hand, generically, only one 



°i.e. disregarding any odd order extrema. 
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of the primary homoclinic points Pq" is bounded away from the fixed point 7(V), 
leading to the creation of large lobes with area of O(A) (see 9 ). On the other hand 
the lobes get stretched and elongated, creating a web of small intersecting regions 
(see 0,111). 

Notice that if KAM tori exist, then the stable and unstable manifolds cannot 
intersect them, hence they are either interior or exterior to R r k 

Lemma 4. (Neishtadt ) Let #„ min denote the phase at which fj,(R(6)) is minimal 
(/i(i?(6*, imin )) = min n(R(0) ) . Then, as uj(v) — > 0, the area enclosed by the largest 
KAM torus inside Rq" (for any n), fJ>(R core (i>)) asymptotes ^{R(dn mil ,))-' 

Furthermore, it is estimated in [H^ that the error scales as 0(ZJ 2 |ln<U| 2 ). 
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